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ABSTRACT 

In most engineering applications, complex states of stress and strain are 
introduced into components during service. It follows that analysis of such 
components requires material descriptions, or constitutive theories, which 
reflect the tensorial nature of stress and strain. In most applications 
involving low homologous temperatures and stress levels below yield, material 
response is linear and elastic. These material characteristics lead to sim- 
ple constitutive relationships involving two material constants. 

For applications involving stress levels above yield, the situation is more 
complex in that material response is both nonlinear and history dependent. 
Plasticity theories have been developed to model these features by using 
flow laws and evolutionary laws expressed in differential form. At elevated 
temperatures, the situation is more complicated still in that material 
response is time dependent. This has led to the development of viscoplastic 
constitutive theories which introduce time by expressing the flow and evolu- 
tionary equation in the form of time derivatives. Most recently, attempts 
have been made to extend viscoplastic theories to transversely isotropic 
materials. The intent here is to develop models which can be used to analyze 
high-temperature components manufactured from advanced composite materials. 

In parallel with the theoretical studies, considerable effort has been 
directed at developing multiaxial testing techniques to verify the various 
theories. This, work faces two technical challenges. The first is identi- 
fying the type of experiment best suited to verifying particular theories. 

In the case of plasticity, yield surface testing techniques were developed 
and used successfully to support the theoretical effort. A similar approach 
has recently been proposed to support development of viscoplastic theories. 

In this case, the time-dependent characteristics of elevated temperature 
response are captured by using inelastic strain rate as a measure of inelas- 
tic state and by determining flow surfaces. It is interesting to note that 
this approach can be extended without difficulty to experimental studies of 
composite materials. 

The second challenge is being able to make stress and strain measurements 
to high levels of precision under biaxial loading conditions. Presently, 
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tension-torsion load cells are fully developed, and so stress measurements 
do not present much difficulty. In contrast, the systems available for meas- 
uring biaxial strains need further development. At room and intermediate 
temperatures, foil strain gage rosettes have been used with some success in 
yield surface studies. Attempts to extend this approach to elevated tempera- 
tures have proved less successful because of the unreliability of high- 
temperature strain gages when used in biaxial stress fields. This has led 
to the development of a number of mechanical extensometers designed specifi- 
cally for use in flow surface determinations at elevated temperature. Again, 
this instrumentation can be used to characterize advanced composite materials 
provided these materials can be fabricated in the form of thin-walled tubes. 

The primary aim of this paper is to outline recent progress in the develop- 
ment of constitutive theories both from the theoretical and experimental 
viewpoints. One important aspect of this work is that it is leading to veri- 
fied material descriptions for advanced composite materials which can be 
implemented in general purpose finite element codes and used for practical 
design. 
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POSTER PRESENTATION 


BIAXIAL EXPERIMENTS SUPPORTING DEVELOPMENT OF PLASTICITY THEORIES 

Yield surface experiments have been used extensively to investigate multiax- 
ial deformation behavior at room and intermediate temperatures. In these 
experiments, tubular specimens are loaded under computer control at fixed 
ratios of torsional stress and axial stress. The specimens are instrumented 
with strain gage rosettes which provide decoupled measures of axial strain 
and tensorial shear strain. During the initial stages of loading, straight 
line relationships are established between stress and strain for both compo- 
nents of loading. The stress levels are increased proportionately until the 
stress-strain response deviates a predetermined amount from linear. These 
deviations, or small offset definitions of yield, are kept small, 25 pe or 
less, to avoid changing the material's state significantly. After achieving 
the target offset, the specimen is unloaded, and the procedure is repeated 
for stress ratios in all quadrants of tension-torsion stress space. The com- 
binations of torsional stress and axial stress giving the required yield con- 
dition are noted at the termination of each probe and used to construct yield 
surfaces . 


INDIVIDUAL LOADING PROBES ARE PROBING SEQUENCE USED TO DETERMINE 

TERMINATED WHEN A TARGET VALUE OF YIELD SURFACES IN TENSION/TORSION 
EQUIVALENT INELASTIC STRAIN F IS STRESS SPACE 
REACHED 


TORSION 
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PLASTICITY THEORIES FOR INITIALLY ISOTROPIC MATERIALS 

The aim of plasticity theories is to provide mathematical descriptions of 
material response at stress levels above yield. Here, material response is 
both nonlinear and history dependent. In the case of initially isotropic 
materials, material descriptions should not be directional. This requirement 
can be met by using as many as three invariants to introduce stress into the 
theory. Classically, the single invariant used for this purpose is the J2* 
or Von Mises, form of equivalent stress. One approach adopted in treating 
history— dependent material response is to use an internal state variable 
ot£j which provides a measure of inelastic state. Evolutionary equations are 
used to relate changes in otij to incremental changes of inelastic strain 
(Prager) or effective deviatoric stress (Ziegler). It is important to note 
that these theories can only be verified through careful experimentation. 


FLOW LAW: dc P, = X d a u 

•I Say 3a k , 

POSSIBLE EVOLUTIONARY EQUATIONS: 


(DRUCKER) 

da,j = C d€jj (PRAGER) 

dofjj = (Sjj — ofjj) d/t (ZIEGLER) 


POSSIBLE YIELD FUNCTIONS: f(J 2 ) = J 2 -K 

f(J2,J 3 ) = J 2 + ^J3 2/3 -K 


(V ON MISES) 
(DRUCKER) 


WHERE 

J 2 = 1/2 SjjSjji J 3 = 1/3 SjjSjkSfcj 

S,j = < 7 jj - 1/3 a kk Sjj; ajj - INTERNAL STATE VARIABLE 
AND c, X, AND n ARE MATERIAL CONSTANTS. 
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COMPARISON OF EXPERIMENTAL YIELD SURFACES WITH THEORETICAL PREDICTIONS 


The yield surface concept described earlier provides a convenient means of 
interpreting plasticity theories. Adopting this approach, the yield function 
describes the shape of the yield surface, and the evolutionary equation 
describes the translation of the yield surface in stress space. Interpreta- 
tion of the data can be simplified further by using modified stress space 
o^2 versus This is because the Von Mises yield function plots as 

a circle in this stress space. Theoretical and experimental yield surfaces 
for a series of nonproportional loadings beyond initial yield are shown here. 
The experiments were conducted at 20 # C on 9Cr-lMo steel and used a 25-ye- 
offset definition of yield (Ellis, 1985). As is typical for most structural 
alloys, the Von Mises yield function gives a close representation of initial 
yield behavior and a poor representation of subsequent yield behavior. Also, 
Prager and Ziegler forms of evolutionary equation provide reasonable approxi- 
mations of yield surface translation in stress space. 


TRANSLATION OF YIELD SURFACES RESULTING FROM NONPROPORTIONAL 
LOADINGS BEYOND INITIAL YIELD 
IDEALIZED EXPERIMENTAL 
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ELEVATED-TEMPERATURE BIAXIAL TESTING SUPPORTING THE DEVELOPMENT 
OF VISCOPLASTIC CONSTITUTIVE THEORIES 


Attempts to extend yield surface testing techniques to elevated temperatures 
have proved less than successful. One major difficulty has been the unrelia- 
bility of high-temperature strain gages when used in biaxial stress fields. 
This led to the development of a number of biaxial extensometers which have 
recently shown promise in probing-type experiments (Ellis, 1983). Another 
difficulty is that the traditional concept of yield breaks down at elevated 
temperatures when material response becomes time dependent. One experimental 
approach developed to resolve this difficulty is to use inelastic strain rate 
eP as a measure of change in inelastic state. By adopting this approach, 
flow surfaces can be determined for particular values of eP and used to 
guide the development of viscoplastic constitutive theories. 


INDIVIDUAL LOADING PROBES ARE 
TERMINATED WHEN A TARGET VALUE OF 
EQUIVALENT INELASTIC STRAIN RATE IP 
IS REACHED 


PROBING SEQUENCE USED TO DETERMINE 
FLOW SURFACES IN TENSIONfTORSION 
STRESS SPACE 

TORSION 


TORSION 



* 12 = ‘12 “ ff 12/2G 
AND ?P = [(e ? 2 ) 2 + 
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A VISCOPLASTIC CONSTITUTIVE THEORY FOR INITIALLY ISOTROPIC MATERIALS 


In addition to treating the complexities of plasticity, viscoplastic consti- 
tutive theories attempt to model the time-dependent features of material 
response at elevated temperature. These features include creep, relaxation, 
recovery, and rate dependence. A convenient method of introducing time into 
the material description is to express the flow and evolutionary equations 
in the form of time derivatives. The majority of viscoplastic theories use 
two internal state variables, <X£j and K, to treat the complexities of 
elevated temperature behavior (Freed, 1988). As in the case of plasticity, 
stress is introduced into the material description by means of invariants. 
Again, the Von Mises form of equivalent stress is used in the majority of 
these models for initially isotropic materials. The validity of this 
approach remains to be verified experimentally. 


FLOW LAW: e[j = 0(T)Z^5f 

EVOLUTIONARY LAWS: ^ = ej - ^ 

H *1 K 2 

j-jVS-OTrW 

WHERE Z(E 2 ), r(K), AND 0(T) ARE MATERIAL FUNCTIONS; H AND h ARE 
MATERIAL PARAMETERS; AND 


Ey = Sjj - 

3 1 

s ij ~ g ff ij “ 2 ° kk 5| i 

r2 1 1/2 

•p [2 .p .pi 1/2 

C 2 = [3 € lj c i| 
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DETERMINATION OF FLOW SURFACES AT ELEVATED TEMPERATURES 


The feasibility of using flow surfaces to investigate multiaxial response at 
elevated temperatures has been demonstrated in a series of experiments con- 
ducted on type-316 stainless steel at 20 and 650 °C (Battiste and Ball, 

1986). As indicated earlier, the main challenge in conducting these experi- 
ments is strain measurement. Obvious requirements in probing-type experi- 
ments are linearity and high resolution. A requirement for minimum crosstalk 
is unique to biaxial testing. This problem arises when loading in the axial 
sense produces apparent torsional strains and vice versa. Further complica- 
tions arise in conducting these experiments at elevated temperatures. The 
problem here is that the electrical noise produced by most heating systems 
causes problems with test system control and resolution. As demonstrated by 
the results shown here, mechanical extensometers have been developed which 
meet most of the requirements discussed previously. Further improvements are 
necessary, however, before some of the more fundamental theoretical assump- 
tions can be tested with a high degree of confidence. 


FLOW SURFACES FOR SMALL VALUES OF EQUIVALENT INELASTIC STRAIN RATE 
CAN BE DETERMINED IN. POST-TEST ANALYSIS OF THE EXPERIMENTAL DATA 


?»■ 


LOADING SEQUENCE 16 
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TIME 
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PRELIMINARY FLOW SURFACES DETERMINED FOR TYPE-316 STAINLESS STEEL AT 20 *C 
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EXTENSION OF BIAXIAL TESTING AND CONTINUUM THEORIES TO 
ADVANCED COMPOSITE MATERIALS 


It has been recognized for some time that the thin-walled tube and tension- 
torsion loading provide an ideal means of investigating deformation behavior 
in composite materials. One advantage of this approach is that properties 
can be determined for compressive loadings without too much difficulty. Fur- 
ther, shear properties can be obtained with minimum ambiguity by using the 
torsional component of loading. Thus, it appears that the high-temperature 
instrumentation and experimental techniques developed for isotropic materials 
can be used to advantage in developing deformation theories for advanced 
composite materials. 


APPLICATION OF CONTINUUM THEORIES TO STRUCTURED MATERIALS REQUIRES 
IDENTIFICATION OF A CONTINUUM ELEMENT (D) THAT IS SMALL COMPARED TO 
CHARACTERISTIC STRUCTURAL DIMENSIONS BUT LARGE COMPARED TO 
CELL SIZE DIMENSIONS 



MEASURES OF ANISOTROPY, u> AND ij, CAN BE OBTAINED BY DETERMINING FLOW 
SURFACES FOR TUBES REINFORCED IN THE LONGITUDINAL AND TRANSVERSE SENSES 
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VISCOPLASTIC CONSTITUTIVE THEORY FOR TRANSVERSELY ISOTROPIC MATERIALS 

Procedures are well established for extending continuum mechanics concepts to 
composite materials. Here, the challenge is to incorporate directionality 
into the material description without violating the laws of mechanics. A 
viscoplastic constitutive theory developed by using this approach is shown 
below for the case of transversely isotropic materials (Robinson et al., 
1987). This form of directionality was incorporated into the theory by spec- 
ifying a strong direction, di. This increased the number of invariants 
initially involved in the theoretical development to seven. This number sub- 
sequently, was reduced to three by tailoring the theory to match assumed 
material response. The strength of the continuum .mechanics approach is that 
it makes use of more than three decades of progress in the field of visco- 
plasticity. This almost assures that the theory will prove successful in 
predicting time dependence and history dependence in composite structures 
involving a single strong direction. Whether the theory can be extended to 
predict the response of complex composite structures is open to question. 

FLOW LAW: «J = f(F)r„ 

EVOLUTIONARY LAW: ajj = h (G) t jj - y (G) 

WHERE 

f(F), h(G), AND 7(G) ARE MATERIAL FUNCTIONS, AND 


F = — j 

K T 2 L 


Ii + 


l ± + ? il_- 

ij | 2 4(4w 2 - 1) 3 . 


r ij = =ij - * [<Wjk + djVia - a o d i d il - 2 fl o (3d i d i " ¥ 


. 1 [_/ I; 

' ~ if 2 r 1 + 

L v 


2 + 2 J 3 

ij 2 4(4w 2 — 1) 3 


CO- 88- 33314 


*i\ = a ij — ^ l d k d i a jk + d j d k a ki 

-2l 0 d i d j ]--fl 0 (3d j 

„ 2 -i 

4(« 2 - 1) 

« = 2 

V 

r_ 4« 2 -1 

V - K l /K t 

« = y l /y t 

Il=J 2 -I+Jl 3 

I 2 =I-I 3 

I 3 =do) 2 

j 2 = 2 S 'i E 'i 

1= d j d jEjkEki 

Io = d j d jEjj 


NOTE THAT THE PRIMED INVARIANTS IN G AND ARE OBTAINED BY REPLACING 
Ejj BY a,j IN THE ABOVE EXPRESSIONS. 
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IMPLEMENTATION OF THE VISCOPLASTIC CONSTITUTIVE THEORY FOR TRANSVERSELY 
ISOTROPIC MATERIALS IN THE MARC FINITE ELEMENT CODE 


The subject viscoplastic model has been implemented in the MARC finite ele- 
ment code and exercised in a number of trial calculations (Arya, 1987). 

The material properties used in these calculations were similar to those 
of a metal matrix composite material, tungsten-copper, at about 800 # F. 
Initially, the analyses were limited to uniaxial stress states. Under these 
conditions, it was shown that the model can treat cyclic plasticity, rate 
dependence, and creep for a range of fiber orientations. More recently, a 
thick-walled cylinder was analyzed to determine performance in calculations 
involving multiaxial stress states. As indicated below, the theory was 
successful in predicting creep-like behavior which was highly dependent on 
fiber orientation. For example, creep in the hoop direction can be seen to 
be negligible for the case of <t> = 90. 


THE SUBJECT CONSTITUTIVE THEORY HAS 
BEEN IMPLEMENTED IN THE MARC FINITE 
ELEMENT CODE AND USED TO ANALYZE A 
THICK-WALLED CYLINDER WITH A RANGE OF 
FIBER ORIENTATIONS (*) 



WHEN LOADED UNDER INTERNAL PRESSURE, 
*4000 psi, THE MATERIAL EXHIBITED 
CREEP-LIKE BEHAVIOR WHICH WAS HIGHLY 
DEPENDENT ON FIBER ORIENTATION 
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